Evidence of Galaxy Cluster Motions with the Kinematic Sunyaev-Zel'dovich Effect 

Nick Hand/'^ Graeme E. Addison,'^ Eric Aubourg,^ Nick Battaglia,^ Elia S. Battistelli,^ Dmitry Bizyaev/ 

J. Richard Bond,^ Howard Brewington/ Jon Brinkmann/ Benjamin R. Brown, ^ Sudeep Das,^°'^^ 

Kyle S. Dawson/^ Mark J. Devlin, ^■^ Joanna Dunkley,^ Rolando Dunner,^^ Daniel J. Eiscnstcin,^'"' 

Joseph W. Fowler, 1*5 Megan B. Gralla,^'^ Amir Hajian,* Mark Halpern,^* Matt Hilton,^^ Adam D. Hincks,*^'^" 

Renee Hlozek,^ John P. Hughes, ^^ Leopoldo Infante, ^'^ Kent D. Irwin, ^^ Arthur Kosowsky,^' ^^ 

Yen-Ting Lin,^^ Elena Malanushenko,^ Viktor Malanushenko,^ Tobias A. Marriage, ^^ Danica Marsden,^'' 

Fehpe Menanteau,^^ Kavilan Moodley,^^ Michael D. Niemack,^^ Michael R. Nolta,^ Daniel Oravetz,^ 

Lyman A. Page,^" Nathalie Palanque-Delabrouille,^^ Kaike Pan,^ Erik D. Reese, ^^ David J. Schlegel,^^ 

Donald P. Schneider, ^^' ^^ Neelima Sehgal,^ Alaina Sheldon,^ Jon Sicvcrs,^"'® Cristobal Sifon,^^ Audrey Simmons,^ 

Stephanie Snedden,^ David N. Spergel,^ Suzanne T. Staggs,^" Daniel S. Swetz,^^ Eric R. Switzer,^ 

^^ Hy Trac,^ Benjamin A. Weaver, ^^ Edward J. Wollack,^° Christophe Yeche,'^® and Caroline ZunckeP^ 

^^ ^Department of Astronomy, University of California, Berkeley, CA 94720 USA 

^SJ ^Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544 USA 

I '^Sub-department of Astrophysics, University of Oxford, 

'—{ Denys Wilkinson Building, Keble Road, Oxford 0X1 3RH UK 

1—5 ^APC, University of Paris Diderot, CNRS/IN2P3, CEA/IRFU, 

t — ^ Observatoire de Paris, Sorbonne Paris Cite, France 

^vj ''Department of Physics, Carnegie Mellon University, Pittsburgh, PA 15213 USA 

^Department of Physics, University of Rome "La Sapienza, " Piazzale Aldo Moro 5, 1-00185 Rome, Italy 

I 1 ''Apache Point Observatory, P.O. Box 59, Sunspot, NM 88349 USA 

^^ ^Canadian Institute for Theoretical Astrophysics, 

r\ University of Toronto, Toronto, ON M5S 3H8 Canada 

. ^Department of Physics and Astronomy, University of Pittsburgh, Pittsburgh, PA 15260 USA 

t£^ "^Department of Physics, University of California, Berkeley CA 94720 USA 

Wh ^'Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720 USA 

^ ^'^ Department of Physics and Astronomy, University of Utah, Salt Lake City, UT 84112 USA 

J, ^'"'Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104 USA 

"*^ ^'^ Departamento de Astronomia y Astrofisica, Pontificia Universidad Catolica de Chile, Santiago 22 Chile 

j^ "'''Harvard College Observatory, 60 Garden Street, MS 20, Cambridge, MA 02138 USA 

i__i ^^NIST Quantum Devices Group, 325 Broadway Mailcode 817.03, Boulder, CO 80305 USA 

^^ Department of Physics and Astronomy, The Johns Hopkins University, Baltimore, MD 21218 USA 

vN '^^ Department of Physics and Astronomy, University of British Columbia, Vancouver, BC V6T 1Z4 Canada 

^ ^'''School of Physics and Astronomy, University of Nottingham, University Park, Nottingham, NG7 2RD UK 

^ ' ^"Joseph Henry Laboratories of Physics, Princeton University, Princeton, NJ 08544 USA 
^^ Department of Physics and Astronomy, Rutgers University, Piscataway, NJ 08854 USA 
^^ Pittsburgh Particle Physics, Astrophysics, and Cosmology Center, University of Pittsburgh, Pittsburgh PA 15260 
'''^Institute of Astronomy and Astrophysics, Academia Sinica, Taipei, Taiwan 

CO '^'^ Department of Physics, University of California, Santa Barbara, CA 93106 USA 

^^ ^^Astrophysics and Cosmology Research Unit, School of Mathematical Sciences, 

Cm University of KwaZulu- Natal, Durban 404^, South Africa 

^—1 ^^CEA, Centre de Saclay, IRFU, 91191 Gif-sur-Yvette, France 

►^ ^^ Department of Astronomy and Astrophysics, The Pennsylvania State University, University Park, PA 16802 USA 

• i-H ^^ Institute for Gravitation and the Cosmos, The Pennsylvania State University, University Park, PA 16802 USA 

j^ '''^Center for Cosmology and Particle Physics, New York University, New York, NY 10003 USA 

^ ''"Code 553/665, NASA/Goddard Space Flight Center, Greenbelt, MD 20771 USA 

<^ (Dated: February 27, 2013) 

Using high-resolution microwave sky maps made by the Atacama Cosmology Telescope, we for 
the first time present strong evidence for motions of galaxy clusters and groups via microwave back- 
ground temperature distortions due to the kinematic Sunyaev-Zel'dovich effect. Galaxy clusters are 
identified by their constituent luminous galaxies observed by the Baryon Oscillation Spectroscopic 
Survey, part of the Sloan Digital Sky Survey Iff. We measure the mean pairwise momentum of 
clusters, with a probability of the signal being due to random errors of 0.002, and the signal is 
consistent with the growth of cosmic structure in the standard model of cosmology. 

PACS numbers: 98.52.Eh, 98.62.Py, 98.70.Vc, 98.80.Es 

Introduction. The growth of cosmic structure over the history of the Universe inevitably results not only in the 
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formation of dense objects, but also in motions of these 
objects. Measurements of these motions have the poten- 
tial to provide both a valuable consistency check on the 
standard cosmological model, and also an independent 
route to constraining cosmological parameters and the 
nature of dark energy. 

In 1972, Sunyaev and Zel'dovich realized that a moving 
galaxy cluster, which is largely composed of hot, ionized 
gas in a dark matter potential well, will induce a small 
brightness temperature shift in the microwave radiation 
passing through it. The shift is proportional to both 
the mass in electrons and the line-of-sight velocity of the 
cluster with respect to the microwave background rest 
frame [3 [5]. This kinematic Sunyaev-Zel'dovich (kSZ) 
effect is distinct from the thermal SZ (tSZ) effect, in 
which scattering from the same hot cluster gas creates 
a spectral distortion (see [3] for a review). In high-mass 
clusters (M ~ 10^^ solar masses), the tSZ signal is typi- 
cally a factor of 20 larger than the kSZ signal; however, 
the two signals are comparable for the low-mass clus- 
ters {M ~ 10^^ Mq) which are far more abundant. (For 
brevity, we refer to any object with mass larger than 
10^'^ Mq as a cluster, even though objects below lO^'* Mq 
are usually referred to as "groups".) The tSZ effect from 
large clusters is now regularly observed in blind surveys 
[1H7], but only upper limits for the kSZ effect from indi- 
vidual galaxy clusters have been achieved to date (SHU] ■ 

In this paper, we present clear statistical evidence of 
the motions of galaxy clusters through their kSZ sig- 
nal in arcminute-resolution microwave maps made with 
the Atacama Cosmology Telescope (ACT) [T^]. Lumi- 
nous galaxies are associated with galaxy clusters [HI [15] , 
and we use the Sloan Digital Sky Survey III (SDSS- 
III) Baryon Oscillation Spectroscopic Survey (BOSS) [TB] 
catalog of these galaxies as galaxy cluster proxies, giv- 
ing the sky location and redshift for thousands of po- 
tential clusters. We then treat the effective microwave 
temperature at 148 GHz measured by ACT in the di- 
rection of the cluster as a noisy estimator of the clus- 
ter's line-of-sight momentum, due to the kSZ tempera- 
ture shift from that cluster. Individual cluster momen- 
tum measurements have a low signal-to-noise ratio, but 
we combine a large number of differential measurements 
to obtain estimates of the mean relative momentum of 
cluster pairs in bins of comoving cluster separation. This 
statistic is insensitive to the tSZ signal, galaxy emission, 
and other sky signals in the ACT data. The conventional 
scenario of structure formation, driven by gravitational 
attraction, predicts that any pair of clusters should have 
a slight tendency to be moving towards each other rather 
than away from each other [T71 [TB] , and we see the ex- 
pected signal in our data at a high statistical significance. 

Survey Data Sets. We make use of two astronomical 
survey data sets. The first is a 148 GHz sky map from 
ACT, a dedicated microwave survey telescope in the At- 
acama Desert of Chile. The map covers a strip approxi- 



mately 3° wide and 110° long with an angular resolution 
of 1.4', centered on the celestial equator and obtained 
over three observing seasons from 2008 to 2010 [13]. We 
match-filter the ACT 148 GHz map with a characteristic 
filter scale at the map resolution of 1.4' [TO], to suppress 
noise from the primary microwave background fluctua- 
tions. Filtered map pixels are 0.5' square, and have a 
noise per pixel ranging from 15 to 25 /iK brightness tem- 
perature; the map is calibrated to 2% by comparing with 
WMAP [20] • Imperfect knowledge of the beam profile 
gives an additional 1% calibration uncertainty on galaxy 
cluster scales. (A similar map at 218 GHz has higher 
noise and is used only for Table 1 below.) 

The second data set is a catalog of luminous galax- 
ies from BOSS Data Release 9 (a combination of the 
CMASS and LOWZ samples from DR9), a component 
of the Sloan Digital Sky Survey HI [^nfej . The cata- 
log contains 27291 galaxies in a 220-square-degree region 
overlapping the ACT sky region (right ascension range 
—43° to -f45°). Galaxies are selected to lie at least 1' 
away from any radio source in the 1.4 GHz FIRST ra- 
dio catalog f24J ; radio contamination is not a significant 
issue. Spectroscopic redshifts range from z = 0.05 to 
z = 0.8 with a mean redshift of 0.51. Luminosities are 
estimated as in Ref. |19j . with an additional mean lumi- 
nosity evolution correction as given in Ref. [5S]. A halo- 
model correlation function analysis shows that most of 
the BOSS galaxies reside in haloes with masses around 
lO^'^ solar masses, with around 10% to 15% in haloes as 
large as lO" M© [26] , 

To estimate the microwave temperature distortion Ti 
associated with galaxy i, we follow the procedure used 
in Ref. [T^]: a 10' by 10' submap centered on the galaxy 
is repixellized into 0.0625' subpixels, convolved with the 
ACT beam profile to smooth the map, and then averaged 
over all subpixels within 1' of the galaxy. The 1' binning 
radius maximizes the signal-noise ratio of our kSZ detec- 
tion, but varying the binning radius between 4" and 4' 
only marginally changes the detection significance. 

The most luminous galaxies in our catalog have the 
largest halo masses. To confirm this, we divide our 
sample into five luminosity bins; Table 1 displays the 
mean central temperature distortion corresponding to 
the galaxies in each bin. The rightmost column gives the 
tSZ distortion brightness temperature at 148 GHz, with 
the next-largest component due to dust emission pro- 
jected out ^27j : this is obtained from a linear combination 
of ACT 148 and 218 GHz signals. Galaxies in the three 
highest luminosity bins, corresponding to about 20% of 
the total, show mean temperature decrements consistent 
with halo-model cluster masses [53] and with the mean 
temperature decrements found in Ref. |19j . 

Mean Pairwise Momentum. Combining the above sur- 
vey data provides a set of galaxy cluster sky positions 
and redshifts (from the luminous galaxy positions and 
redshifts) and line-of-sight momenta (from the ACT tem- 
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TABLE I: Mean brightness temperature fluctuations in sky directions corresponding to the BOSS DR9 galaxies, in bins of 
galaxy luminosity. The right column corresponds to the tSZ brightness temperature at 148 GHz, STtsz = '5ri4g — 0.325 (5T2i8, 
projecting out a dust emission component |27|. The ACT maps are match-filtered at an angular scale of 1.4', equal to the beam 
size at 148 GHz, then subpixelized, convolved with the beam profile, and summed over all subpixels within 4" of the galaxy. 



perature). To compare this data set with cosmological 
models, consider the mean pairwise momentum statistic: 



PpairCr) = ((Pj 
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where galaxy cluster i has momentum p; and comoving 
position Vi, the comoving separation vector between a 
pair of clusters i and j is r^j = r, — r^, overhats denote 
unit vectors, and the average on the right side of the 
equation is over all cluster pairs in a bin around comoving 
separation r = \rij\. If two galaxy clusters are moving 
towards each other, their contribution to Ppair('') will be 
negative, and if moving apart, positive. An estimator of 
Ppair(^) using only line-of-sight momenta is [28] 
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where 9 is the angular separation between two clusters on 
the sky and r^ = |ri| is the comoving distance to cluster 
i, which can be computed from the cluster redshift us- 
ing standard ACDM cosmological parameters [53]. (The 
cluster velocity gives negligible contribution to the dis- 
tance estimate for clusters at a cosmological distance, and 
a bias in estimating Eq. (pi) which is small compared to 
the measurement errors.) The statistic Ppa,ir{r) is equal to 
the familiar mean pairwise velocity ■ypair(^) [50^5^ times 
the average mass of the clusters in the sample. 

We can measure the line-of-sight component of the mo- 
mentum via the kSZ microwave temperature fluctuation, 
Tksz.j = — ^ksz Prfjj assuming that the ratio of the total 
cluster mass to its mass in hot gas is simply the univer- 
sal ratio of matter density Q^ to baryon density fib |33| . 
The normalization A^ksz depends on the pixel scale and 
beam size of the microwave map, and the cluster density 
profile. Simulations including these effects [31j give an 
expected mean temperature signal in the ACT 148 GHz 
map of 1.6 and 0.3 /iK for clusters with a typical line- 



of-sight velocity of 200 km/s and masses 10^'* and 10^^ 

The statistic Ppair('") is both linear and differential, 
giving it desirable systematic error properties [35j . Any 
microwave temperature signal associated with individual 
galaxy clusters, like the tSZ effect, will average to zero 
as long as it does not depend on the relative distance 
between cluster pairs. Redshift-dependent signals can 
contribute to Ppair(^) f^nd be confused with the cluster 
kSZ signal, including infrared emission from galaxies in 
the cluster (which increases with redshift out to z = 2), 
any radio source emission, and small variations of the 
tSZ signal due to evolution of average cluster mass and 
temperature. However, we can measure these effects on 
average by simply finding the average microwave tem- 
perature T{z) corresponding to clusters at a given red- 
shift, and correcting the temperature in the direction of 
an individual galaxy cluster for this redshift-dependent 
piece. We evaluate a smoothed T{z) by averaging the 
temperature towards all galaxies, each with redshift Zi 
and a Gaussian weight factor exp[— (z — Zi)^/2tT^] with 
az = 0.01; our results are nearly insensitive to the value 
of az within a wide range. The resulting T(z) has a mean 
near zero and an absolute value of up to 3/.tK. 

We thus evaluate the mean pairwise kSZ signal, cor- 
recting for possible redshift-dependent temperature con- 
tributions, as 
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This quantity differs from Eq. (pi) by the amplitude fac- 
tor A^ksz- Figure [1] displays this statistic for the ACT 
pixel temperatures corresponding to the 5000 most lumi- 
nous BOSS DR9 galaxies in the ACT sky region {L > 
8.1 X 10^° L0); this luminosity cut minimizes the total 
noise from combined Poisson and pixel noise. Also dis- 
played is the signal extracted from a kSZ-only sky sim- 
ulation, based on underlying large-volume cosmological 
simulations |34j . adjusting the mass limit of the simu- 
lation halos to give the best fit to the data. We infer 



that our galaxy luminosity cut corresponds to a cluster 
halo mass limit of roughly M200 — 4.1 x 10^^ Mq and a 
mean cluster halo mass of M200 = 6.5 x 10^'' Mq- Error 
bars are estimated via bootstrap resampling. Neighbor- 
ing bins have a mean correlation of 0.25 and we include 
smaller mean correlations out to a 5-bin separation, as 
determined using independent simulation volumes. 

The measured points largely fall below zero and have 
Ax^ — 23 for 15 degrees of freedom, compared to the 
best-fit model. The model is a good fit to the data: 13% 
of random data realizations with the same normal er- 
rors and correlations have larger A^^- The measured 
points have A^^ = 43 for 15 degrees of freedom, com- 
pared to a null signal; the probability of random noise 
having Ax^ at least this large is 2.0 x 10^'^ including 
correlations. The measured points approach zero signal 
as the comoving pair separation increases, which demon- 
strates that the signal depends on spatial separation, not 
redshift separation. 

Null tests are simple, as the statistic is essentially a 
sum of pixel temperatures, half with positive and half 
with negative signs, with weights corresponding to rel- 
ative galaxy positions. Figure [T] also displays the null 
test corresponding to using the same weights but random 
pair positions compared to the signal plot (A^^ = 11.6 
for 15 degrees of freedom). Success of this null test veri- 
fies that the function T{z) correctly models any redshift- 
dependent contributions to the microwave signal. Chang- 
ing the sign in the second term of Eq. Q from negative to 
positive also gives a null signal (Ax^ = 9.9 for 15 degrees 
of freedom) . 

Discussion and Prospects. The signal in Fig. [T] repre- 
sents the first measurement of the cosmic velocity field 
made directly with respect to the rest frame of the Uni- 
verse. It is consistent with simulations based on the stan- 
dard cosmological model. This signal is also the first 
clear evidence for the kinematic Sunyaev-Zel'dovich ef- 
fect. A recent attempt by Kashlinsky et al. to measure 
the large-scale bulk flow via the galaxy cluster kSZ sig- 
nal uses galaxy clusters from X-ray surveys and searches 
for an overall dipole dependence of the microwave tem- 
perature in the WMAP data at these locations [36l |37] . 
However, Keisler [35] found the first reported detection 
was not statistically significant. Osborne et al. [39] rean- 
alyzed the most recent results including both a monopole 
and dipole term, obtaining limits on a bulk flow a factor 
of three below the reported detection of Ref. [37| . Mody 
and Hajian [ID] also fail to reproduce the bulk flow re- 
sult using Planck and ROSAT galaxy clusters. Planck 
will soon make a more precise test of this reported large- 
scale flow [H]. The statistic used in this paper is differ- 
ential, which mitigates many of the potential systematic 
errors affecting bulk flow measurements, but also is not 
sensitive to an overall bulk flow. 

Most previous work on peculiar velocities using opti- 
cal observations has measured the properties of the local 
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FIG. 1; The upper panel shows the mean pairwise momentum 
estimator, Eq. Q, for the 5000 most luminous BOSS DR9 
galaxies within the ACT sky region (red points), with boot- 
strap errors. The solid line is derived from numerical kSZ sim- 
ulations [33] using a halo mass cutoff of M200 =4.1xlO^'^M( 
The probability of the data given a null signal is 2.0 x 10 
including bin covariances. The lower panel displays the same 
sum but with randomized map positions, and is consistent 
with a null signal. 



bulk flow, but has not been able to extend measurements 
to cosmologically interesting distances. The traditional 
method of measuring velocities - a Doppler shift of an 
object's radiation spectrum - is very challenging at cos- 
mological distances because the spectrum of an object is 
redshifted due to the expansion of the Universe, and this 
cosmological redshift is typically large compared to the 
velocity frequency shift. Precise distance measurements 
are required, a difficult observational problem. Recent 
optical work [12] extends to around 100 Mpc, a redshift 
of 2; = 0.02, while this paper uses galaxy cluster velocities 
out to z = 0.8. Future large optical surveys such as the 
LSST may enable competitive cosmological velocity sur- 
veys using large catalogs of standard candles for distance 
measurements [33]. 

The evidence for a nonzero mean pairwise momentum 
from a kSZ signal presented here can also be interpreted 
as a measure of baryons on cluster length scales; a deficit 
of observed baryons has long been a cosmological puz- 
zle [45. Our signal is roughly consistent with the stan- 
dard baryon fraction based on primordial nucleosynthe- 
sis, given independent halo mass estimates based on clus- 
tering of our luminous galaxy sample. This issue will be 
addressed in a future paper. 



Future improved measurements of the mean pairwise 
velocity have the potential to put strong constraints on 
dark energy and modified gravity [45ll47j . The measure- 
ment we have presented here is the first step on a new 
path to constraining structure growth in the Universe. 
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